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Abstract 

The hadronic decays rj — > tttttt, rj — > 7T7r7r and rj — > r/7T7r are investigated within 
a £7(3) chiral unitary approach. Final state interactions are included by deriving 
the effective s-wave potentials for meson meson scattering from the chiral effective 
Lagrangian and iterating them in a Bethe-Salpeter equation. With only a small set 
of parameters we are able to explain both rates and spectral shapes of these decays. 
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1 Introduction 



The hadronic decays r] — > tttttt and rf — > 7T7T7r are of great interest since they violate isospin 
symmetry. At such low energies the calculations of the decay amplitudes are based on 
effective chiral Lagrangians which have the same symmetries and symmetry breaking 
patterns as the underlying QCD Lagrangian. They are written in terms of effective 
degrees of freedom, usually the octet of pseudoscalar mesons - pions, kaons and the eta - 
[T], however, they can easily be extended to include the rf also [2J. With the effective 
Lagrangian at hand, one can perform a perturbative expansion of the decay amplitudes in 
powers of the pseudoscalar meson masses and external momenta. For the decay 77 — > tttttt 
the expansion parameters are the ratios of the Goldstone boson octet masses or Lorentz 
invariant combinations of external momenta over the scale of spontaneous chiral symmetry 
breaking, m|/A^ and (pp -pp/)/A 2 , respectively, with A x = 4:irf n ~ 1.2 GeV and higher 
loops contribute to higher chiral orders. At low energies these ratios are considerably 
smaller than unity and the chiral expansion is expected to converge. 

The inclusion of the rf, on the other hand, spoils the conventional chiral counting 
scheme, since its mass does not vanish in the chiral limit so that higher loops will still 
contribute to lower chiral orders. This can be prevented by imposing large N c counting 
rules within the effective theory in addition to the chiral counting scheme. In the large 
N c limit the axial anomaly vanishes and the rf converts into a Goldstone boson. The 
properties of the theory may then be analyzed in a triple expansion in powers of small 
momenta, quark masses and 1/N C , see e.g. [2JIH1IHE!- In particular, m v > is treated as a 
small quantity. Phenomenologically, this is not the case since = 958 MeV, and we 
will therefore treat the rf as a massive state leading to a situation analogous to chiral 
perturbation theory (ChPT) with baryons. Recently, a new regularization scheme - the 
so-called infrared regularization - has been proposed which maintains Lorentz and chiral 
invariance explicitly at all stages of the calculation while providing a systematic counting 
scheme for the evaluation of the chiral loops Infrared regularization can be applied 
in the presence of any massive state and has been employed in chiral perturbation theory 
including the rf (U(3) ChPT) in [7||HJin], but its applicability is restricted to very small 
external three momenta. This is surely not the case for the decay rf — > tttttt where large 
unitarity corrections are expected due to final state interactions between the three pions. 

Final state interactions have already been shown to be substantial in 77 — > tttttt both in 
a complete one-loop calculation in SU (3) ChPT JUj and under the inclusion of final state 
interactions using extended Khuri-Treiman equations [TT]. The last calculation comes 
closest to the experimental decay rate but still remains below it, if the usual quark mass 
ratios are employed. However, the experimental rate of the decay can be reproduced 
by increasing the quark mass ratio [(rrid — m u )/{m s — m)] ■ [{rrid + m u )/{m s + rh)] from 
1/(24. 1) 2 to 1/(22. 4±0.9) 2 which is now considered to be the most accurate determination 
of this quark mass ratio. It is therefore obligatory to include final state interactions also 
in the decay rf — > tttttt. We will apply the following approach in the present investigation: 
after deriving the effective potentials for meson meson scattering from the chiral effective 
Lagrangian, we iterate them utilizing a Bethe-Salpeter equation (BSE). This method 
has been proven to be useful both in the purely mesonic sector and under the inclusion 
of baryons ^21 The BSE generates dynamically bound states of the mesons and 
baryons and accounts for the exchange of resonances without including them explicitly. 
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The usefulness of this approach lies in the fact that from a small set of parameters a large 
variety of data can be explained. It has been extended recently to U(3) ChPT in [T4"| 
where effects of the rf in meson meson scattering are investigated and the parameters 
of the U(3) chiral effective Lagrangian are constrained by comparing the results with 
the experimental phase shifts. Again, with only a few chiral parameters the phase shifts 
have been reproduced. Using the same technique, the electroproduction of the rf meson 
on nucleons has also been investigated yielding good agreement with data ^Hj- This 
method allows us to account for final state interactions and to study the importance 
of resonances in the different isospin channels, a topic of interest, e.g., in the dominant 
hadronic decay mode of the rf , rf — > tjittt. A full one- loop calculation of this decay 
has been performed in [Oj utilizing infrared regularization and reasonable agreement with 
experimental data was achieved. However, both higher order final state interactions and 
the explicit inclusion of resonances have been omitted. It has been claimed, on the 
other hand, that this decay can be described in a tree-level model via the exchange of 
the scalar mesons <r(560), / (980) and a (980) which are combined together with k(900) 
into a nonet [Ej (see also [13 EH E3 E01 El])- The authors find that the exchange of 
the scalar resonance ao(980) dominates. By iterating the effective chiral potentials to 
infinte order in a Bethe-Salpeter equation, we can investigate the importance of these 
resonances explicitly. Hence, the present work provides a unified approach to the decays 
7] — > 7T717T, rf — > tttttt and rf — > tjtttt with chiral symmetry and unitarity being the main 
ingredients. With only a few chiral parameters we can compare our results with a variety 
of experimental data and even predict the decay rate for rf — ► 7r°7r + 7r~. For simplicity, we 
will restrict ourselves to s-waves since good agreement with experimental data is already 
obtained; an extension to higher multipoles is straightforward. 

In the next section, we present the chiral effective Lagrangian up to fourth chiral order 
and the ifir)rf mixing for different up- and down-quark masses m u ^ ma is discussed in 
detail as it arises from the second and fourth order Lagrangian, respectively. In Sectional 
the implementation of the final state interactions via the BSE is illustrated. In Sections 0] 
to 14.11 the decays r\ — > nnn, rf — > txtxtx and rf — > rjirn are discussed. Our results are 
compared with experimental data. 

2 The Lagrangian density 

In this section the Lagrangians at second and fourth chiral order are presented and the 
resulting 7r°-?7 8 -?7o mixing is discussed. The effective Lagrangian for the pseudoscalar 
meson nonet (ir, K, rjg, rjo) reads up to second order in the derivative expansion Ej 

where U is a unitary 3x3 matrix containing the Goldstone boson octet (ff , tt°, K ,K°, rjs) 
and the i] . Its dependence on tt°,7] 8 and r/o is given by 

U = exp(diag(l, -1, 0) ■ in°/F + diag(l, 1, -2) ■ it] 8 /V3F + iV2t] /V3F + ...). (2) 

The expression (...) denotes the trace in flavor space, F is the pion decay constant in the 
chiral limit and the quark mass matrix M. = diag(m M , m<f, m s ) enters in the combination 
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X = IBM. with B = — (0\qq\0)/F 2 being the order parameter of the spontaneous sym- 
metry violation. As we do not consider external (axial-) vector currents in the present 
investigation, the covariant derivatives have been replaced by partial ones. 

The coefficients Vi are functions of 770, Vi(f] /F), and can be expanded in terms of this 
variable. At a given order of derivatives of the meson fields U and insertions of the quark 
mass matrix M one obtains an infinite string of increasing powers of 770 with couplings 
which are not fixed by chiral symmetry. 4 Parity conservation implies that the Vi are all 
even functions of r/ except V 3 , which is odd, and 14(0) = V 2 (0) = 14(0) - 31^(0) = \F 2 
gives the correct normalizaton for the quadratic terms of the mesons. The potentials Vi 
are expanded in the singlet field r] 



V 

13 



= v? ) +vW$ + vP$ + ... for , = 0,1,2,4 

*i .4 1) | + 4 3) | + ... 0) 



with expansion coefficients to be determined phenomenologically. 

In order to describe the isospin violating decays 77 — ► tttitt and rj' — > mm, we need to 
distinguish between the up- and down-quark masses, m u and m d , which leads to TT°-rjs-r)o 
mixing. Taking the Lagrangian from Eq. (fl|). diagonalization of the mass matrix to first 
order in isospin breaking m u — m d yields the mass eigenstates 7t ,t},t] / with 

7f° = 7T° - eq - 2edrf 
r/s = e7r° + r] — 'drj 

Vo = 3e^vr° + $7} + 7]'. (4) 

The angles e and d arise from the mass differences of the light quark masses m u , ma, m s . 
While breaks SU (3) symmetry due to different strange and nonstrange quark masses, 
e is proportional to strong isospin violation — m u 

V3m d -m u m 2 f 



m s -m a/3 ( 



ml — ml 



^^ B(ro ,- A )^ K -„4) (5) 



where we used the abbreviations 

«?> = \F 2 - \^vf 
m\ = B(m d - mu) 

m=\(m d + m u ). (6) 

The quantity m\ can be expressed in terms of physical meson masses by applying Dashen's 
theorem [22], which implies the identity of the pion and kaon electromagnetic mass shifts 
up to 0(e 2 p 2 ) 

2 2 2,2 2 frr\ 

m e = m R0 - m K± + m n± - m n0 . (7) 
4 Note that we do not make use of 1/N C counting rules. 
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At leading order in the chiral expansion the meson masses are given by 

ml = 2Bm, 
m 2 K ± = B(m s + m u ), 

m 2 K0 = B(m s + m d ). (8) 

Isospin violation is known to be small, hence e is a small quantity despite being of zeroth 
order in the quark masses. Terms of order (m^ — m u ) 2 are tiny and therefore neglected 
throughout. In the isospin limit of equal up- and down-quark masses m u = m^ the angle 
e vanishes and 7f° does not undergo mixing with r] 8 and 770 • 

However, this is not the whole story. As shown in [E] the fourth order Lagrangian 
contributes to 77-7/ mixing at leading order as well which is due to the fact that we count 
the 7/ mass as quantity of zeroth chiral order. The fourth order Lagrangian has the form 

£ (4) = Ei A (770)0* (9) 

with the fourth order operators 

4 =-(C^)(M), 5 = -(C»C»M), 

17 = (C»)(C,){M), 18 = -{C»){C,M), 

6 =(M)(M), 7 =(N)(N), [W) 

8 =l(MM + NN), Oi2 = \{MM - NN), 

We used the abbreviations = Wd^U, M = U ] x + X ] U, N = U ] x ~ X ] U and the 
coefficients can be expanded in 7/0 in the same manner as the Vi in (JHJ). The mixing 
between fc°, rj 8 and 770 cannot be described in terms of just two angles any longer. Due to 
the operators in Eq. ()10|) two subsequent transformations of the original fields vr°, 77 8 and 
770 are necessary to bring both the kinetic and mass terms of the Lagrangian into diagonal 
form. At leading order in isospin breaking the new transformed fields are then related to 
the original ones by 



7T 



1 + R^o)tx q + R^^rj + Rjto v 'V' 



f] 8 = R 87T o7c° + (1 + R 8v )v + R 8 - n 'f]' 

r]o = Ro n on° + R Qv r] + (1 + Rq v ')v' (H) 



with the mixing parameters given by 

p(0) _ m l p(0) _ _ p(0) 

V3(m^ - ml) 



R 8w° ~ ^iw°(^iLo + 3^Gmo), R jr°r)~ R 8t!° i^Sr, + 3 A 



GMO J 



(2) _ ^\m 2 q - ml) (2) _ (2) 8^ 0) 18 (m 2 v - m 2 ] 



(2) _ op(0) p(2) p(2) _ 9 p(0)p(2) 

(°)™2 , «/o(0)/ 2 , _2N nfl(O) / 2 



(12) 



R (2) CX + 6^K + "4) p( 2) 2^>^ 18 (m 2 +mi 

K ir ( >TT» ~ p2 ' ""'Or/' ~~ ^2 



(2) 4^m|+6^(m2+m 2 ' 

^8r/ — ^ 2 
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where the superscript denotes the chiral order and we have employed the abbreviations 
0$Ui 8 = WP+ffi-WW+WFl and Pi% = (3^+3(3^/2. The quantity m 2 = 2v^ / ' F 2 
is the mass of the 77' meson in the chiral limit and the deviation from the Gell-Mann- 
Okubo mass relation for the pseudoscalar mesons is given by 



Am 2 K - ml - 3m 2 Q(m 2 - 



2\ r a/~k 

(°)_Rfl(0)_1Ofl(°) 



'GMO 9 9 t-,9 



^-6^-12^+ /l2;/;2 



(13) 



Finally, the remaining terms of the fourth order Lagrangian that we need in this 
calculation are given by 

O =(CC v C ll C u ), 1 =(C»C li )(C v C v ), 

2 =(C*C V )(C IX C V ), 3 ={C»C li CC v ), 

13 = -(C»)(C,C»C U ), C> 14 = -(C»)(C,)(C»C U ), 

15 = -(C»)(C»)(C,C U ), 16 = (C»)(C,)(C»)(C U ), (14) 

2 i = (CC^iN), 22 = (CCJiiN), 

O23 = (CHCpiN), 24 = {C»)(Cr)(iN), 

25 = (zMN), 26 = (M)(iN). 

Usually the /3q term is not presented in conventional ChPT, since there is a Cayley- 
Hamilton matrix identity that enables one to remove this term leading to modified co- 
efficients Pi, i = 1,2,3,13,14,15,16 [I], but for our purposes it turns out to be more 
convenient to include it. Hence, we do not make use of the Cayley-Hamilton identity 
and keep all couplings in order to present the most general expressions in terms of these 
parameters. One can then drop one of the /3j involved in the Cayley-Hamilton identity at 
any stage of the calculation. 



2.1 Values of LECs 

The unknown coupling constants of the chiral Lagrangian - the so-called low-energy con- 
stants (LECs) - need to be fit to experimental data. This has already been accomplished 
in where after applying the same approach as in the present investigation we con- 
strained the LECs of the Lagrangian up to fourth chiral order by comparing the results 
with the experimental phase shifts of meson meson scattering. Agreement was achieved 
in the isospin I = 0, | channels up to 1.2 GeV and in the isopin I = |,2 channels up 
to 1.5 GeV. The discrepancey to the data above 1.2 GeV in the 1 = channel, e.g., is 
due to the omission of the 4n channel, see j23j- In order to obtain reasonable agreement 
with the experimental phase shifts, the values of the involved chiral parameters are then 
constrained to the following ranges 

/?£ 0) = (0.6 ± 0.1) x 10~ 3 , 4° ) = (-0.5 ± 0.1) x 10" 3 , , , 

pf ] = (1.4 ± 0.2) x 10- 3 , /?< 0) = (0.2 ± 0.2) x 10- 3 . ' :J ' 

while keeping = = \F 2 and all remaining parameters set to zero. However, 
variations of some of the parameters which have been set to zero do not yield any sig- 
nificant effect for the phase shifts, and in principle they can have finite values. Such a 
parameter is (5^ or the combination from Eq. (jSJ), and we have chosen (3^ = and 
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7T° 




Figure 1: Shown is a possible contribution to final state interactions in the decay 
i] — > 7r°7r + 7r-. Here, it takes place in the s-channel between tt + tt~ and the bubble chain 
represents part of the solution to the BSE at energy ^/s. 

= = \F 2 in |T3] but we are free to change their values in the present investi- 
gation as long as they do not contradict the data for the phase shifts. As a matter of 
fact, this work provides a test whether the same values for the LECs can be employed 
when describing the hadronic decays of the 77 and 7/ and will allow for finetuning of the 
parameters that were not constrained in ^3], e.g. f3^ and A good fit to the decays 
discussed in this work is given by 

«?> = = 0, 

/5j 0) = 0.56 x 10- 3 , (3f } = -0.3 x 10" 3 , (16) 
/^ 0) = 1.4 x 10" 3 , pf ] = 0.06 x 10- 3 . 

with t>2 = jF 2 — y/Ev^ — 3v^ and all the remaining parameters being zero. Note 
that there have been moderate changes with respect to the values in Eq. (JT5j) which were 
necessary to obtain better agreement with the experimental decay rates and the spectral 
shapes of the Dalitz distributions. It is surprising that with a small number of parameters 
we are able to explain a variety of data within the approach. 

It is also important to note, that the parameter choice in Eq. (fTtjj) is not unique, since 
variations in one of the parameters may be compensated by the other ones. Nevertheless, 
we prefer to work with this choice; it is capable of describing the experimental phase 
shifts and the hadronic decays of the 77 and 7/ - as we will see shortly - with a minimal 
set of parameters and motivated by the assumption that most of the OZI- violating and 
1 / iV c -suppressed parameters are not important, although (3^ and v% , i n ^2 3^2 
have small but non-vanishing values. 

3 Final state interactions 

An investigation of the hadronic decay 77 — > imn at leading order in ChPT [23 
yields a decay width which is significantly below the measured value j2H]- The result is 
improved by performing a next-to-leading order calculation, but still fails to reproduce the 
phenomenological value. The substantial increase from leading to next-to-leading order 
already indicates that higher order effects beyond one-loop ChPT may be important as 
well and should not be neglected. Kambor, Wiesendanger and Wyler [TT] were able to 
approximate these higher order effects by using Khuri-Treiman equations, which describe 



7 







Figure 2: Solution to the Bethe-Salpeter equation T = A + AG A + AG AG A + ... 

final state interactions. The underlying idea is that the initial particle, i.e. the r\ or rf, 
decays via chirally constrained vertices derived from the effective Lagrangian into three 
mesons and that two out of these three mesons rescatter an arbitrary number of times, 
see Fig. ^for illustration. There are three possible ways in combining two of the mesons 
to a pair while leaving the third one unaffected which corresponds to the s-, t- and it- 
channel, respectively. Interactions of the third meson with the pair of rescattering mesons 
are neglected. Such an infinite meson meson rescattering can alternatively be generated 
by application of the Bethe-Salpeter equation. To this end, the s-wave potentials for 
meson meson scattering are derived from the chiral effective Lagrangian and iterated in 
a Bethe-Salpeter equation. The BSE then generates the propagator for two interacting 
particles in a covariant fashion. 

The Bethe-Salpeter equation for the two-particle propagator T from a local interaction 
A is given bv [TH 127] 

rpf \ a( \ I f iddk T(p,qi,k)A(p,k,q f ) 

T( P , ft , Qf ) = A( P , ft , qf ) + J — d (fc2 _ m2)(fc2 _^ 2) - (17) 

Combinatorial factors for two identical particles in the loop have to be taken into account, 
but we prefer to keep this form of the BSE and modify the amplitudes accordingly. In 
that case A must be the four-point amplitude from ChPT multiplied by — ~ in order for 
T to be proportional to a bubble chain with the correct factors from perturbation theory. 
The factor of \ is the symmetry factor of two identical particle multiplets in a loop and 
— 1 stems from factors of i in the vertices and propagators. 

The solution for T amounts to an infinite resummation of the interaction kernels A 
and the free propagators G and can be diagrammaticaly depicted as an infinite interaction 
chain, Fig. 121 In a short notation, both the equation and its solution can be written as 

T = A + TGA, T= {1-GA)- 1 A = A + AGA + AGAGA + ... (18) 

For the decays considered here we would like to use these final-state interactions in such 
a way as to approximate to a large extent the one-loop result from conventional ChPT. 
The second term in the expansion of T, AGA, equals the unitary corrections in one- 
loop ChPT (i.e. without tadpoles). Tadpoles which are not included in our approach 
have been shown to yield numerically small effects in scattering processes in the physical 
region and can furthermore be partially absorbed by redefining chiral parameters |12j . 
We have therefore neglected the tadpoles and crossed diagrams and restrict ourselves to 
the tree diagrams from next-to-leading order ChPT in the calculation of the potential A. 
As mentioned above, for the decays one must take into account the three possible ways 
in combining two of the light mesons to a pair while leaving the third one unaffected, 
which corresponds to the s-, t- and u-channel. Final state interactions will occur in all 
of these channels and in order to obtain the full amplitude one must add an interaction 
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chain in each channel, T s + T t + T u . This reproduces the unitary corrections at one-loop 
order; the leading contribution in this sum, however, does not yield the tree level result 
A from ChPT. Portions of the tree level amplitude are included in each T and must be 
subtracted from the bubble chain. The proper expression can be easily accounted for by 
choosing the full amplitude as T = A + (T s — A s ) + (T t — A t ) + (T u — A u ) which includes 
the unitary corrections in one-loop ChPT, while the leading contact term pieces A Si t, u are 
removed from the solutions T Syt)U and replaced by the tree diagram amplitude A. 

In the evaluation of the BSE we will restrict ourselves to s-waves, since they are 
expected to dominate in the decays discussed here, and indeed this simplification gives 
results already in agreement with experiment. An improvement of the results - particu- 
larly for the spectral shape of the Dalitz distributions - may be expected by including the 
p-waves, but this is beyond the scope of the present investigation. We can further simplify 
the integral in the BSE (|17|). since we are only interested in the physically relevant piece 
of the solution T with all momenta put on the mass shell. The amplitude A contains in 
general off-shell parts which deliver via the integral a contribution even to the on-shell 
part of the solution T. However, these off-shell parts yield exclusively chiral logarithms 
which - besides being numerically small - can be absorbed by redefining the regulariza- 
tion scale of the loop integral. Furthermore the off-shell parts are not uniquely defined 
in ChPT. We will therefore set all the momenta in the amplitudes in ()17|) on-shell, which 
has also been done in previous work, see e.g. Finally, we need to consider all two- 

meson channels. In particular, those which have the same set of quantum numbers can 
couple amongst each other. Consequently, the objects A, T, G are promoted to matrices 
in channel space and are functions of the two particle invariant mass squared p 2 . With 
these simplifications the BSE equation reduces to (fTSj) as a matrix equation for every p 2 
and the solution is obtained by matrix inversion. The matrix G is a diagonal matrix with 
elements given by the integral 

mm[P ] = J (2^p (k 2 -m 2 ){{k-p) 2 -m 2 ) (19) 

where m and m are the masses of the two mesons in the corresponding channel. In 
dimensional regularization the integral is evaluated as 



16tt 2 



, mm m 2 — m 2 , m 
1 + In — =- + - In — 

H z p A m 



^\l^rafh{p 2 ) t h \/^rnfh{p 2 ) 



p 2 (m + m) 2 — p 2 

Xrnrnip 2 ) = {{m - mf - p 2 ) ((m + m) 2 - p 2 ) . (20) 

The integral is divergent and \i takes the role of a regularization constant which can be 
chosen independently for every channel. With these simplifications we write the full four- 
point amplitude T a b c d as a sum of the solutions T a b tC d of the BSE in Eq. (fTTj) in the s-, t- 
and li-channel 

T a bcd(s, t, u) = A abcd (s, t, u) + (T- A) ab>cd (s) + (T - A) acM (t) + (T — A) adM (u), (21) 
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where, e.g., (T — A} a & jCC j(s) is a bubble chain in the s-channel with particles ab on one 
side, particles cd on the other side, and the tree level piece removed. The expression 
(T — A) ab cd (s) is thus equivalent to the s-wave projection of the final-state interactions 
of T, whereas the tree-level amplitude A abcd (s,t,u) is not decomposed into partial waves. 

In the considered decays, rj, rj — > tttttt and rj — > r]7T7r, the two-particle states are either 
uncharged or simply charged. There are nine uncharged combinations of mesons 

ttV, tt+tT , t]tt°, TjTj, K° K° , K + K~ , r/vr , 7/77, 7/7/ (22) 

which are labeled with indices 1, . . . , 9, respectively, and a set four charged channels 

TT°TT + ,r]TT + ,K + K°,r]'TT + (23) 

labeled with indices 1,...,4. Due to isospin breaking, couplings between the different 
channels of a set occur, but charge conservation prevents transitions between both sets. 
One obtains a 9 x 9 matrix A which summarizes the amplitudes for the uncharged 
channels and a 4 x 4 matrix A + for the charged channels. The division of the two-particle 
states into charged and uncharged channels simplifies the treatment of the BSE. We have 
presented the leading order contributions for A and A + in the appendix. 

In order to obtain the final expressions T abc d for the decay amplitudes, the solutions 
T° and T + of the BSE in Eq. (jl7j) must be corrected due to the overall symmetry factor 
—1/2 which was initially absorbed into the matrix A. This is achieved by setting 

7^0,00 = —y/%T® x , T^QfiQ = —\[2T® 1 , 

T TO T TO 

J »yO,H — — — J 32' J r;'0,H — — — J 72' 

T v +,o+ — —T 2 i, T v i +>0+ = —T 41 , (24) 

Tri>ri,oo = ~V2Tg 1 , TrfQ^o = — Tg 3 , 

rp rpQ rp rp + 

1 V' r l^ — — J 82 5 1 r/' + ,V+ ~ -M2' 

where an index 0, ± denotes ir ,^. The same relations hold between the amplitudes 
A ab)Cd and Ap°. The additional factors a/2 in some of the channels are included in order 
to account for the statistical factor occuring in states with identical particles 



4 T] —> 7T7Y7Y 

The decay r] — > tttttt violates isospin and can take only place due to a finite mass differ- 
ence m u — m d or electromagnetic interactions. The latter ones are expected to be small 
(Sutherland's theorem) [2H1 which has been confirmed in an effective Lagrangian frame- 
work [TI]. Disregarding these, the amplitude is proportional to m u — m d and provides a 
sensitive probe on SU (2) symmetry breaking. 

The process rj — > tttttt is the dominant decay mode of the rj with the measured rates 
being j2Hj 

T(r] -> 7r 7r 7T ) = (379 ± 40) eV 

r(?7 -> tt°tt + tt-) = (274±33)eV (25) 



10 



and the ratio given by 



T 77 -> 7T 7T 7T , , 

r = — M n - ; = 1.35 ±0.05 26 

T(r] -> 7T°7T + 7T-) V 7 

Further information is provided by the energy dependence of the Dalitz plot which 
must be compared with our results. To this end, we introduce a compact notation which 
can be applied to all decays considered in the present work. Let <j>\ — > 020304 be the 
decay of a meson with mass M into three lighter mesons with masses m 2 , 777,3 an d w,±. For 
the cases considered here, we will choose 777 3 = 7774 = 777.34. The Dalitz plot distribution is 
conventionally described in terms of the two variables 

V3(u - 1) 

x - 



2M(M - 777 2 - 2m 34 ) 



(m 2 + 2m 3i )[(M-m 2 ) 2 -s] _ j 
V 2Mm u (M - m 2 - 2m 34 ) ^ ' 

with the Mandelstam variables 

8 = (pi- P2) 2 , t = (px- P3) 2 , w = (pi - Pi) 2 (28) 

with pk and rrik being the four-momentum and mass of the k-th meson and p\ = M 2 . The 
Dalitz plot distribution \A(x, y)\ 2 is then expanded as 

\A(x,y)\ 2 = N[l + ay + by 2 + cx 2 ]. (29) 

for the charged decay and 

\A(x,y)\ 2 = N'[l + g(y 2 + x 2 )] (30) 

for the neutral decay with no linear term in y and only one coefficient g due to Bose 
symmetry. The expansion parameters have been determined experimentally to be 



Layter et al. [2H]: a = -1.08 ± 0.014, b = 0.034 ± 0.027, c = 0.046 ± 0.031 

Gormley et al. jSOJ: a = -1.17 ± 0.02, b = 0.21 ± 0.03, c = 0.06 ± 0.04 

Amsler et al. [SI]: a = -0.94 ± 0.15, b = 0.11 ±0.27, 

Amsler et al. [22]: a = -1.19 ± 0.07, b = 0.19 ± 0.11, 

Abele et al. [33]: a = -1.22 ± 0.07, b = 0.22 ± 0.11, c = 0.06 fixed 

Aide et al. 133: g = -0.044 ± 0.046, 

Abele et al. [35]: g = -0.104 ± 0.04, 

Tippens et al. [3H]: g = -0.062 ± 0.010. 

(31) 

4.1 Tree diagram contributions to 77 — > tttttt 

The first step in our approach is to calculate the tree diagram contributions with vertices 
from the Lagrangians in Eqs. (0E]). The tree diagram contributions for A v _> n o n + n - from 
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£(°+ 2 ) and £< 4 ) read 



3\/3F2(m2 - to 2 ) 



x 8m 2 (F^-F^) 2A gmo (8to 2 - 3s) 



3(to 2 — rn^)F^ 3(4to 2 — 3s) 

(32) 



)fft + u - s) 2 - s 2 - ft - u) 2 ) 

F 2 (4m 2 -3s) ll + J V u ) ) 



where we have employed the expressions for the meson masses and the decay constants 
at fourth chiral order 8J. It is interesting to see that this expression from U(3) ChPT 
equals the one obtained in ^Hj in SU(3) ChPT. The differences due to rjrj' mixing are 
completely absorbed by a change in Aq M o. The expression for the neutral decay follows 
immediately via the AI = 1 selection rule 

Ajooo 0, t,u) = A,o+- (s, t, u) + A70+- (t,u,s) + A,o+- (u,s,t). (33) 

Taking the central values for the LECs from Eq. (|T6*j) and applying Dashen's theorem 

B(m d - m u ) =m 2 e = m 2 K o - m 2 K+ + m 2 + - m 2 , (34) 

we obtain for the rates of the decays 

T(rj -> tt°ti + tt-) = 88 eV, T(r] -> ttVtt ) = 134 eV (35) 

with the ratio given by r = 1.52. Clearly, the results are far from being even close to 
the experimental value; this is well-known for quite some time |24[ l2*K] . A full one-loop 
calculation within conventional ChPT PU], although improving the result considerably 
(r = 167 ± 50 eV), still fails in being in agreement with the current data. The unitarity 
corrections, particularly in the tttc channel, turn out to be large and indicate that a 
summation of two-pion rescatterings is necessary which cannot be done in a loopwise 
expansion. Higher order corrections to Dashen's theorem as discussed in the literature [SZ1 
EH EH] may increase the result a bit further, but cannot account for the large discrepancy 
to the data. We will return to this point in the next section. 
The resulting parameters of the Dalitz plot distribution are 

a = -1.14, 6 = 0.31, c = 0.014; g = 0. (36) 

While the parameters for the charged decay seem to be in agreement, the coefficient g 
of the uncharged decay vanishes identically in this approximation and is in disagreement 
with the recent precise measurement of 



4.2 Final state interactions in 77 — > 717171 

As already mentioned in the introduction, final state interactions due to two-pion rescat- 
terings in the decay 77 — > nun have been evaluated in 11, using Khuri-Treiman equations 
and making use of the one-loop ChPT result, see also |4T)| I4*T] for a similar approach. 
Our framework is much simpler than that and provides a unified aproach to the hadronic 
decays of the rj and rj'. The underlying idea is that the initial particle, in the present case 
the rj, decays via chirally constrained vertices derived from the effective Lagrangian into 
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three mesons and that two out of these three mesons rescatter. There are three possible 
ways in combining two of the mesons to a pair while leaving the third one unaffected 
which corresponds to the s-, t- and w-channel, respectively. 

Following (|2*TJl and (j24j) in our approach, the full amplitudes read 

T„ 0+ -(M,«) = Afl+-(8,t,u) - (T - A)° 32 (s) — (T — A)+(t) — (T — A)+(u) (37) 



T v0O0 (s,t,u) = A v0Q0 (s,t,u) - V2(T - A)° 31 (s) - V2(T - A)° 31 (t) - V2(T - A) 



3l(«) 



Note that the AI = 1 selection rule, Eq. ()33)1 . is not valid any longer, since our method 
includes higher orders in the isospin breaking parameter m u — which spoil this relation. 
The difference to this rule is therefore proportional to higher orders in m u — m d , and we 
verified numerically that the deviations are rather small. Using the parameters from 
Eq. (fH)j) we obtain the results 

T(r] -> 7i\ + 7i-) = 237eV, T(r] -> ttW) = 351 eV (38) 

with the ratio given by r = 1.48. Our decay rates are close to the experimental ones, and 
for the neutral decay we are even within the error bars. 

We have employed Dashen's theorem Eq. ([34)1. in order to arrive at the results Eq. ([38)1 . 
however, slightly modified values for the quark mass ratios increase the result even further 
[TT] . To this end, one replaces the quark mass difference m u — m d by 

B(m d - m u ) = 7L— 7T (m 2 K - m 2 n ) (39) 

with 

q2 = m s -m m s +m 
m d - m u m d + m u 

Dashen's theorem yields the value 

2 2 2 

^Dashcn = ^ ^ _ m jl + + m 2 + _ ^ = ^ ' ( 41 ) 

Alternatively, the decay rj — > 3n can be used to determine Q. In order to reproduce the 
experimental values of Eq. (|25j) we obtain 

Q = 23.4 ±0.8, (42) 

a value which is in agreement with the requirement by Dashen's theorem and indicates 
that higher order corrections to this low-energy theorem may be small, in contradistinction 
to what has been discussed in [^7[ EE1 EH] ■ 

Let us turn now to the comparison of the spectral shapes for these decays. Our method 
yields 

a = -1.24, 6 = 0.37, c = 0.05; # = -0.014. (43) 

The results are in the same ballpark as those from jTT] and indicate that our method 
captures the same physics as the elaborate work of ^T]. The situation is more complex 
when our results are compared with the data in Eq. (j31|) . The parameters a and b are in 
agreement with |321 , [HH] (a) and (b), respectively, but not with the other ones. Also, 
our value for g is excluded by the most recent measurement [HE] > but since the experimental 
uncertainties are significant in most cases, a new determination of the spectral shape with 
higher statistics is needed. 
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4.3 Resonances in 77 — > tttttt 

It is interesting to investigate which resonances could be of importance for the decay and 
yield a substantial contribution to the decay width. We are able to study the effects of 
the resonances within our framework, since we can in principle switch on and off the final 
state interactions (j2~T|) in the s, t, u channels separatly. For the neutral decay, however, 
this would destroy Bose symmetry and we will restrict ourselves to the charged decays, 
for which we will switch on the t- and w-channel only simultaneously, in order to preserve 



Bose symmetry. The results 


for the decay 77 




tx 7x + tx - 


are 








tree: 


T = 88 eV, 


a 


= -1.14, 


b = 


0.31, 


c = 


0.014 


tree + s-channel: 


T = 316 eV, 


a 


= -1.04, 


b = 


0.23, 


c = 


0.007 


tree + 1, it-channel: 


T = 50 eV, 


a 


= -1.53, 


b = 


0.61, 


c = 


0.093 


tree + s, t, it-channel: 


T = 237 eV, 


a 


= -1.24, 


b = 


0.37, 


c = 


0.005 



(44) 



Two of the three decay products can form either an isoscalar or an isovector state and their 
invariant mass must lie between Im^ pa 275 MeV and m v — m n pa 410 MeV. Furthermore, 
the t, it-channels are charged, hence the 1 = resonances cannot propagate in these 
channels. The so-called a resonance is in the appropriate mass range, while the resonances 
/o(980) and a (980) may also influence the decay, although their masses are comparatively 
large. From (jUj) one observes that final state interactions in the t, li-channel decrease the 
decay width by a small amount. In the resonance picture this could be interpreted as the 
tail of the 1 = 1 resonance ao(980). Large effects on the decay width stem from s-channel 
final state interactions where also 1 = resonances are allowed. The /o(980) is as far 
away from the physical region as the a (980) and is thus expected to have a small impact 
on the decay, so that the enhancement in the decay width can be mainly attributed to 
the a. 



5 rf —> TTTTTT 

We will now apply the same method to the decay rj' — > txtxtx in which the final state 
interactions should also play an important role. According to our knowledge, there exist 
only tree level calculations for this decay in which unitarity corrections have been ne- 
glected throughout [121 E31 El EEJ • All of them disagree with the present experimental 
knowledge. E.g., in the framework of large N c ChPT this decay was calculated up to 
next-to-leading order [33] (loops start at next-to-next-to-leading order) yielding the de- 
cay rates T(rj' — > 7T 7r 7r°) pa 2200 eV and T(r]' — > 7r°7r + 7r~) pa 1600 eV which are almost 
an order of magnitude larger than the data. The measured rates of this decay mode are 
given by j2El 

r(r/ -> ttVtt ) = (311 ± 77) eV, r(r/ -> ttVtt-) < 1005 eV (45) 
with the ratio 



IX IX IX 1 

>0.3. (46) 



r(l/ — > 1X°TX + 1X ) 

Hence, the theoretical situation is unsatisfactory requiring the inclusion of unitary cor- 



rections. 



14 



5.1 Tree diagram contributions to rf — > unit 

The tree diagram contributions with vertices from the Lagrangians 

£(0+2) j n Eqg> (JJJ ^ 

read for A v / 0+ _ 



4\/6m £ 2 (?B^ — 2ml 



A^ +-(s,t,w) = — — — T . (47) 

9(m^ — m^)F 4 

The pertinent tree level amplitude for the neutral decay follows from the AI = 1 selection 
rule. Employing Dashen's theorem, Eq. (|34j). leads to 

T(rj' -> tiVtT) = 114 eV, r(r/ -> ttVtt ) = 55 eV (48) 

with the ratio given by r = 0.48. As expected, we cannot reproduce the experimental 
rate for the charged decay, but it is worthwile mentioning that our result for the charged 
decay remains below the experimental decay rate in contradistinction to the rates in 
H2 E3 El HH] ranging from 450 eV in jSj to 1600 eV in gSj. 

The Dalitz plot parameters read in the tree diagram approximation 

a = -7.72, 6 = 5.78, c = 0.87; # = 0.94. (49) 

Unfortunaly, experimental data on the spectral shape is not available. 



5.2 Final state interactions in rj — > tttttt 

The smallness of our results in Eqs. (|4*K|) hints at the importance of unitary corrections. 
Applying the same model as in the previous section we arrive at 

T(rj' -> n°n + TT-) = 294eV, T(rf -> 7r 7r 7r ) = 291 eV (50) 

which is in good agreement with the available data. In fact, the result for the charged 
decay is a prediction and could serve as a check of our approach in a future experiment. 
The ratio of the two decays is given by r = 0.99. The slope parameters of the Dalitz plot 
are given by 

a = -3.08, 6 = 0.13, c = 0.62; p = -0.86, (51) 

and need to be verified experimentally. However, these predictions must be treated with 
care, since for these decays we expect contributions from p-waves as well due to the 
presence of the p(770). 



6 fj — > T]7T7T 

The last decay we are going to discuss in this work is the dominant decay mode of the rj', 
i]' — > ?77r7r. The measured rates are |26j : 

r(r/ -> t]tt + tt-) = (90±8)keV, T(rj f -»• j^ttV) = (42 ± 4) keV (52) 

with the ratio given by 

r(r/ — > f]7T + 7T~ 
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2.1 ±0.4. (53) 



Although isospin violating corrections seem to be small for this decay, we will not mod- 
ify our approach and keep different up- and down-quark masses. The experimentally 
measured slope parameters are 

Kalbfieisch a = -0.16 ± 0.06, 

Aide et al. jUj: a = -0.116 ± 0.026, 6 = a 2 /4±0.13 2 , c = 0.00 ± 0.03 (54) 

Briere et al. jH]: a = -0.042 ± 0.050 5 , 



6.1 Tree diagram contributions to r{ — > ijtttt 
The tree level amplitudes are given by 

A,' v +-(s,t,u) = A,„,ooO,t,w) = ^pf (55) 

They lead to the decay rates and slope parameters 

T(rj' ^rin°n ) = 6.9 keV, a = -0.204, b = -0.005, c = -0.060; , . 

r(V ->• r]n + n-) = 12.2 keV, a = -0.194, b = -0.005, c = -0.055. 

Clearly, we are not able to explain the data just by taking into account the tree level 
contribtutions and need to include unitary corrections as for the other decays. 



6.2 Final state interactions in rf — > tjtttt 

By iterating the effective s-wave potentials from meson meson scattering we arrive at the 
decay rates 

r(r/ -> ??7r 7r ) = 43.8 keV, a = -0.105, b = -0.065, c = -0.004; , , 
r(V -> 777T+7T-) = 77.7keV, a = -0.093, 6 = -0.059, c = -0.003 '° 

being in qualitative agreement with experiment and with a relatively fixed ratio of 1.8. 



6.3 Resonances in 7/ — > j]ttti 

We can turn on final state interactions as described in Sec. 14.31 in both the s and t, u 
channels independently to investigate the importance of resonances. The results for the 
neutral decay read 



tree: V 


= 6.9 keV, 


a = 


-0.204, 


b = 


-0.005, 


c = 


-0.060 


tree + s-channel: V 


= 10.0 keV, 


a = 


-0.455, 


b = 


-0.064, 


c = 


-0.050 


tree + t, w-channel: T 


= 37.8 keV, 


a = 


+0.021, 


b = 


-0.005, 


c = 


+0.008 


tree + s, t, w-channel: T 


= 43.8 keV, 


a = 


-0.105, 


b = 


-0.065, 


c = 


-0.004 



(58) 

and for the charged decay 



tree: 


r = 


12.2 keV, 


a = 


-0.194, 


b = 


-0.005, 


c = 


-0.055 


tree + s-channel: 


r = 


17.9 keV, 


a = 


-0.424, 


b = 


-0.059, 


c = 


-0.043 


tree + t, u- channel: 


r = 


66.2 keV, 


a = 


+0.021, 


b = 


-0.004, 


c = 


+0.006 


tree + s, t, w-channel: 


r = 


77.7 keV, 


a = 


-0.093, 


b = 


-0.059, 


c = 


-0.003 



(59) 
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The main contribution to this decay arises from the final state interactions in the t, u- 
channels, i.e. when the rj' decays into a pion and two mesons which then rescatter infinitely 
many times to finally pass into a pion and an eta. The iteration of the effective potentials 
in these channels mimics the effect of the ao(980) which is usually generated within chiral 
unitary approaches, see • We therefore confirm the importance of the a® for this decay 
as it has been claimed, e.g., in jTHj. While in the tree level approach of ^E] the ao has 
been included as an explicit degree of freedom, it automatically results in our approach 
by considering the final state interactions. At first sight, these findings may seem to be 
contradictory to the conclusions in pj, where the importance of the coupling constant 
combination (5^ + (3^ was stressed which dominates this decay via a contact interaction. 
However, in the framework of 9\ effects of the ao are hidden in /?< 0) + (3^ via resonance 
saturation j3H| and contribute indirectly through the contact interaction of (3^ + f3^ . 



7 Conclusions 

In the present work, we studied the hadronic decays rj — > 7T7T7r, rf — > unit and rj' — > rjinr. 
To this end, the most general chiral effective Lagrangian which describes the interactions 
between members of the lowest lying pseudoscalar nonet (ir, K, rj, rj') was presented up 
to fourth chiral order. The pertinent tree diagram amplitudes derived from this effective 
Lagrangian are calculated. In all cases, the results for the decay rates are in sharp 
disagreement with the measured values and suggest that unitary corrections are dominant 
in these decays. 

In order to improve the theoretical situation, we included final state interactions. They 
were modelled by deriving the effective s-wave potentials for meson meson scattering from 
the chiral effective Lagrangian and iterating them in a Bethe-Salpeter equation. The 
underlying idea is that the initial particle, rj or rj', decays via chirally constrained vertices 
derived from the effective Lagrangian into three mesons and that two out of these three 
mesons rescatter. There are three possible ways in combining two of the mesons to a 
pair while leaving the third one unaffected which corresponds to the s-, t- and w-channel, 
respectively. The iteration of the effective potentials in these channels mimics the effect 
of resonances and may be used to confirm calculations in which the same resonances have 
been put in by hand. With only a small set of parameters we are able to explain both 
rates and spectral shapes of these decays. 

The decay rj — > tttttt violates isospin and omitting the small contributions from elec- 
tromagnetic interactions it can only take place due to a finite mass difference m u — m^. 
The amplitude is then proportional to m u — rrid and provides a sensitive probe on SU (2) 
symmetry breaking. Therefore, the n°T]Tf mixing for different up- and down-quark masses 
m u ytz ma is discussed in detail as it arises from the second and fourth order Lagrangian, 
respectively. Our results are in reasonable agreement with experiment after employing 
Dashen's theorem. Conversely, this decay may be used to constrain the double quark mass 
ratio Q and we obtain Q = 23.4 ± 0.8, a value which is in agreement with the require- 
ment by Dashen's theorem and indicates that higher order corrections to this low-energy 
theorem may be small. 

For the decay rf — > unit which also violates isospin we obtain agreement with the 
available data. Our result for the neutral mode, T(rf — > 7r 7r°7r°) = 291 eV, is within the 
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experimental error ranges T(rj' — > 7r°7r°7r ) = (311 ± 77) eV, while we can even deliver 
a prediction for the charged mode, T(r]' — > 7r°7r + 7r~) = 294 eV, for which only an upper 
experimental limit of 1005 eV exists. The slope parameters for the spectral shapes of the 
Dalitz plots may also serve as a check of our method. 

The third and final decay we discussed is r/ — > 7/7T7T. Again, reasonable agreement with 
data was achieved and we could furthermore confirm within our approach the importance 
of the ao (980) resonance for this decay as has been claimed in the literature. This is 
also consistent with the findings in [9 , where effects of the ao are hidden via resonance 
saturation in the contact interaction of the coupling constant combination (3^ + /3g which 
dominates the decay. 

We were able to show that the results are in overall good agreement with the decay 
rates, although the approach includes only a few parameters. In the case of the slope 
parameters we are also in the right ballpark, but the experimental situation is less settled 
here and one needs new experiments with higher statistics in order to reduce the error 
bars. Since we included only s-waves in our approach, we cannot expect, of course, to 
reproduce the spectral shapes of the decays very precisely. The p-waves, although small 
in magnitude, will alter our results for the Dalitz parameters, and need to be included in 
order to make more accurate statements. 

The method introduced in the present work has also set the stage for investigation 
on further 7] and rj' decays. Among these are the anomalous decays 77/77' — > 77 and 
vM ~ > 7r + 7r~~7, the decays 77/77' — > 7r°77, and the strong CP violating decays rj/i]' — > 7r7r. 

A Amplitudes 

In the appendix we present the bare s-wave amplitudes A that are iterated in the Bethe- 
Salpeter equation (|T5j) and related to the four-point functions at the tree level by equations 
analogous to (|24jl . We present here for brevity only the contributions from the leading 
order Lagrangian £( 2 \ Contributions from C^> are treated in a similar way. 

A.l Uncharged channels 

The amplitude matrix A for the nine uncharged channels consists of an isospin conserving 
part A° f0 and an isospin breaking part i.e. A = A° e0 + A° el . The matrices A° e0 and A° el 
can be further decomposed as 



The submatrices B and D describe the couplings between the channels (7r°7r°, 7r + 7r~,777r°, 
r]r],K K , K + K~) and {rfn , T]'r], rft]'), respectively. Transitions between the first chan- 
nels and the latter ones are given by the matrix C and its transpose. The matrices B, C 



and D read (with vf ] = \F 2 - fvW^ - ^ + \^vf ] and = \F 2 - 2y/bo$ ) + 




1 = 



0,1. 



(60) 
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A. 2 Charged channels 

For the four charged channels (ir°ir + , r]Tc + , K + K°, r]'Tc + ), the full matrix A + is given by 
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